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ABSTRACT 



The standard unbroken electroweak theory is known to erase baryon number. 

The baryon number symmetry can be restored in the core of a neutron star 

as its density diverges via gravitational instability due to a binary merger 
^ . event. We argue that for certain double Higgs models with discrete symmetries, 

OC . this process may result in an expanding self-sustained burning front which 

would convert the entire neutron matter into radiation. This process would 
^ . release ~ 10^"^ ergs in electromagnetic radiation over ~ 10~^ sec, with negligible 

r^ . baryonic contamination. The resulting fireball would have all the properties 

--».^ , necessary to produce a 7-ray burst as a result of its interaction with ambient 

p ! interstellar gas. The subsequent Higgs decay would produce a millisecond burst 

of ~ 10^^ ergs in ~ 100 GeV neutrinos which might be observable. The above 



O 

d '. giving rise to the observed matter-antimatter asymmetry today. 



mechanism may have also caused electroweak baryogenesis in the early universe. 



H ' Subject headings: cosmology: theory - gamma rays: bursts 



Introduction 



The energy source for 7-ray bursts (ORB) has been enigmatic since their discovery 
(see Fishman & Meegan 1995, for a recent review). The recent detection of absorption 
lines at a redshift z = 0.835 in the optical afterglow spectrum of the event GRB970508 
(Metzger et al. 1997), established a firm lower bound on its distance, and hence on its total 
radiation energy ^ 10^^ erg (Waxman 1997). The long duration of some GRBs is difficult 
to reconcile with their millisecond time structure. If the variability is caused by an external 
shock due to the interaction of the ffieball with an ambient medium (Meszaros & Rees 
1993a,b), then the conversion efficiency must be small (Sari & Piran 1997a) and the total 
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energy released approaches the rest mass energy of a neutron star. A further challenge to 
GRB models is the baryon contamination problem. The high Lorentz factors implied by the 
short variability and non-thermal spectra of GRBs (see, e.g. Woods & Loeb 1995) require 
that the baryonic fraction be ^ 1%, of the total energy. This constraint is particularly 
difficult to satisfy in stellar environments such as neutron star mergers (Ruffert et al. 1997, 
and references therein), where a considerable amount of baryonic debris is unavoidable. 

The standard electroweak theory allows for baryon number violating processes ('t 
Hooft 1976). Through a combination of chirality and topology, baryon and lepton number 
may be violated. A change in the winding number of the gauge fields results in fermions 
being pulled out or pushed into the Dirac sea. Normally, in the broken symmetry vacuum, 
such effects are very strongly suppressed since the gauge fields are massive and difficult to 
excite. However, in the unbroken phase it is well established that when the electroweak 
symmetry is restored, the baryon and lepton number violation processes are unsuppressed 
(for a review, see Rubakov & Shaposhnikov 1996). It is natural to ask whether these 
intriguing processes could be excited in the Universe today. There are two fundamental 
barriers which make their appearance rare: (i) a large energy input is required to change 
the gauge/Higgs winding numbers, and (ii) the energy must be input in a coherent fashion 
into the long wavelength components of the bosonic fields (and this, for example, does not 
happen in high energy collisions of single elementary particles). 

In this Letter we argue that baryon number symmetry might be restored in the core of 
a neutron star if its density is raised to extreme values. Such a process might be triggered 
during the catastrophic gravitational collapse caused by some cataclysmic event (such as 
binary coalescence or substantial accretion), which raises the neutron star mass above 
its stability limit. Once the symmetry is restored in a sufficiently large (macroscopic) 
volume, a self-sustained burning front would propagate throughout the star and convert 
the entire neutron matter into radiation, releasing ~ 1 GeV per baryon in the process. 
The burning front will be accompanied by an electroweak domain wall in which the barrier 
to baryon number violation is greatly reduced. The existence of dense nuclear matter 
allows long wavelength modes of the gauge fields to be naturally excited. Overall, the 
burning front could release ~ 10^"^ ergs in electromagnetic radiation with negligible baryonic 
contamination. The resulting fireball would have all the properties necessary to produce a 
GRB due to its interaction with ambient interstellar gas. 

In §^ we present the hydrodynamic solution for the associated electroweak burning 
front. We then discuss the trigger mechanism for the explosive process described above, 
in §^. §§ considers the energetics of electroweak baryon number violation on moving 
walls. Finally, we examine the implications of our model for cosmological baryogenesis and 
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electroweak theory in §1^, and summarize its testable predictions for GRBs in 



2. Self-Similar Explosion 

In our model, we consider three regions of the dynamical explosion. In the outer 
region, the cold neutron star material is at rest. The shock front propagates through this 
material at some speed A^c, where c is the speed of light. This burning front is the second 
region under consideration, which has a characteristic thickness of an electroweak length 
~ 10~^^ cm. The burning front is sustained by the presence of an electroweak domain wall. 
While the microphysics in this region is complicated, we can describe its hydrodynamic 
implications as if it were infinitely thin. We will return to its discussion in ^. The global 
conservation of energy and momentum thus dictate the global dynamics of this burning 
region. The last region under consideration includes the postshock (burnt) radiation phase. 
This heated region will initially expand due to the increase in pressure. This pressure is 
directly exerted against the domain wall, which then propagates further outward. The goal 
of this section is to quantify the self-sustained dynamics of these regions, assuming that the 
process has already been triggered. The possible trigger mechanism will be considered in 
the next section. 

The equations of motion of a relativistic fluid are given by the energy- momentum 
conservation equations T.'f^' = where T'^'^ = (p + p)u^u^ +pg^^ . We use a metric signature 
convention of (— , +, +, +). In the rest frame of the star, we have 

dp + pP"^ 1 d r^{p + p)l3 ^ 



dt 1-/3^ r^dr 1 - (3"^ 
d {p + p)(3 Id r^{p + p)f3^ dp ^ 
dt 1-p^ r^dr l-(3^ dr ' ^ ^ 

where r is the radial coordinate and the speed of light is set to unity. We ignore gravity 
for now, assuming that the star is initially in static equilibrium, and that the final phase is 
gravitationally fully unbound. 

We assume that the pre-shocked neutron star material is described by an equation of 
state po = apo- The pre-shock density of a neutron star, po ~ 2 x 10^'' g/cm^, is slightly 
above nuclear density. The pressure is typically small, a ~ 0.1, but can reach relativistic 
values near the center of the star, where a ~ 1/3. The burnt phase consists of radiation 
and electroweak false vacuum energies, pbumt = Pf + Pv, Pbumt = Pf/3 — pv- Typically 
Pv ~ Po/100, as discussed further in §^. The burning front has a width of order the 
electroweak scale, ~ 10~^^ cm, much smaller than the macroscopic scale of the neutron 
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star, and can therefore be approximated as a sharp discontinuity. Since the neutron 
star density is nearly constant (Shapiro & Teukolsky 1983), the explosion should be well 
approximated by a self-similar solution (see examples in Blandford & McKee 1976; Pen 
1994). The amount of energy released per baryon equals to its rest mass and therefore 
the Lorentz factor of the burning front cannot be much greater than unity. This implies 
that any self-similar self-sustained burning front (which does not decay as it expands) must 
propagate at a constant speed through the neutron matter. 

We therefore seek solutions to (P and (|l]) which depend on r, t only through the 
combination A = r/t. We define the dimensionless densities p = p/po and pv = V{(f)~)/po, 
where V is defined in §^ below. After some algebra, we find that the equation for (3{X) 
separates from the equation for p: 

A/?' = 2/3(1 -/?^)(1-/JA) 

^ /?2(3-A2)-4/5A + 3A2-l' ^ ^ 

and the density is given as a simple integral over the fluid velocity, 

~ ~ i f'^ 8/3 (A- - /3) dX' \ 

^ = ^^^"P1A /?2(3-A-)-4/?A' + 3A--l|- ^3) 

This self-similar scaling imposes a constant shock speed with velocity A^. The subscript 
{...)s on a fluid variable denotes its value immediately behind the shock front, e.g. 
P. = hm;,^;,-p(A). 

In order to solve the system of equations (H) and (||), we need to specify the postshock 
boundary conditions. Equations (|l|) and (|l|) provide the postshock velocity in the star rest 
frame and the postshock density in the fluid frame: 



. _ 2[A,^(1 -pv)-a- Pv] - V-3(l + afXs^ + 4[a - A,'(-l + py) + py] 

^' X,i3- a- Apy) ' ^^> 

^ ^ 3(l + c.)(l-/32)A, 

^' A{Xs-(3s){l-PsXsy ^' 



Requiring the postshock velocity to be real in the shock frame yields the inequality 

3 + a{3a - 2) + 8py{py + a - 1) + v^(l + a)J{Apy + a- 3)(4py + 3a - 1) 
TT^ — < A, < 1 

(6) 
which simplifies to 3/4 < A^ < 1 if a = p^ = in the idealized case of a low-pressure 
neutron star and a negligible burning potential. We conclude that the burning front must 
propagate supersonically in the star frame. In the shock frame, the post-shock velocity 
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of the fluid is smaller than its pre-shock value, therefore this solution corresponds to a 
detonation front. By examining equation (Q) we find that all solutions must pass through 
the point P = 0, \ = l/^/S. The only physical solution in the range < A < l/^/S is /9 = 0, 
so any solution is matched continuously at the boundary with discontinuous derivatives. 
Normally one specifies a boundary condition away from the shock front, but in our case 
this boundary condition is automatically satisfied. The shock speed A^ is therefore a free 
parameter which can lie anywhere in the allowed range. The actual value of As is either set 
by the trigger mechanism (i.e., the initial conditions) or by the viscous forces acting on the 
domain wall boundary. In the generic case where viscous forces slow down the shock to its 
minimum velocity, the unique solution is 

2 _ 3 a pv , 

-^^ " 4 " 2 " T' ^^' 

to leading order in a and pv- The postshock density is n^ ^ 3 — 4q; — 7pv and the exact 
postshock velocity 



,2_^ 2(l + a) 



3 — a — Apv ' 

giving the constraint 3a; + Apv < 1. Even though the velocity derivative 13' formally 
diverges, the solution is well defined. 

By virtue of equation (|1|), all shock solutions automatically satisfy total energy 
conservation 

which provides a consistency check for the final numerical solution. Figure |^ shows the 
numerical solutions for the velocity and density profiles behind the shock for three different 
choices of A^, assuming a = and py = Q. 



3. Trigger Mechanism 

We have found that if triggered, a self-sustained burning front would propagate at 
a constant speed and spontaneously convert baryons at nuclear density into radiation. 
However, we still need to demonstrate how this process might be spontaneously triggered, 
and under which conditions the burning front would release sufficient energy to sustain 
itself. 

The rate of neutron star mergers in the Hubble volume is similar to the observed rate of 
GRBs (Narayan et al. 1991; Phinney 1991). It is therefore natural to examine the question 
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whether an electroweak burning front would be triggered by such mergers. Typically, the 
combined mass of the two neutron stars would exceed the maximum mass for a stable star 
at nuclear densities, and the merger product would be gravitationally unstable and tend to 
collapse to a black hole. An electroweak burning front might be produced if there is an inner 
region where the density is sufficiently high to restore electroweak symmetry before being 
engulfed by a black hole horizon. In this region, baryon-number will be strongly violated 
and most baryons will decay into anti-leptons. Similar conditions might occur in binary 
systems where a neutron star exceeds the maximum stable mass due to accretion from its 
companion. However, the formation of black holes by much more massive systems, such 
as the direct collapse of massive stars, are characterized by much lower densities (at the 
Schwarzschild radius, p oc 1/M^) and are not expected to trigger the electroweak transition 
before their horizon forms. 

The details of the trigger process depend on the equation of state of neutron star 
matter above nuclear density, which is not yet fully understood. Some nuclear equations 
of state predict that p > p/3 for densities about ten times nuclear (Prakash et al. 1997). 
The sound speed cj. = dp/ dp then obtains a maximum value above that of a pure radiation 
fluid. At yet higher densities, the sound speed eventually declines back to the 1/3 value of 
a locally interacting radiation field due to asymptotic freedom of QCD. If we now consider 
a neutron star which collapses slowly through phases of hydrostatic equilibrium, then the 
central equation of state will stiffen until it reaches its maximum sound speed, and then 
soften again. The region which enters the softened phase will collapse faster than its 
surroundings, potentially leading to a runaway density. The compression of this region is 
caused by the weight of the column of matter above it. The electroweak symmetry would 
be restored in this region at a density ~ 10^^ g/cm^, some twelve orders of magnitude 
above the initial density. This region will not be surrounded by a black hole horizon if its 
radius is ^ 1 cm, or its mass is ^ 4 x 10~^Mq. A lower bound on the radius of the trigger 
region arises from the requirement that a self-sustained burning front needs to release more 
energy interior to it than the total energy stored in the domain wall. At nuclear density 
the energy released is ~ (200MeV)^i?^ while the wall energy is ~ (250GeV)^/2^, implying a 
radius ^ 10~^ cm. Higher background densities will only weaken this constraint. 

We have two apparently contradictory requirements. The fluid must collapse to 
electroweak densities, and then expand again to drive the explosion. In the dense phase, the 
relative pressure gain from the burning of each nucleon is small. One would expect that the 
gravitational collapse and the subsequent phase transitions will raise the entropy of the fluid 
in the central region and cause entropy inversion, therefore resulting in convective flows. A 
fluid element taken through the electroweak phase transition which would subsequently rise 
due to buoyancy forces, could trigger the explosion as it cools adiabatically and is trapped 
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in the false vacuum. 

4. Kinematics of Baryon Number Violation on Propagating Walls 

In this section we discuss the energetics and timescales of the baryon number violating 
processes on bubble walls moving through cold neutron matter. We consider the possibility 
that neutrons are converted into antineutrinos as the bubble wall passes. Our analysis will 
be extremely schematic, at best a pointer to the key issues involved. At the end of this 
section we shall comment on the additional work which is needed in order to settle the 
microphysics of electroweak bubble walls in neutron matter. 

We shall consider particle physics models which include more than one Higgs doublet 
field. These models are amongst the most plausible extensions of the standard model and 
occur as sub-models of low energy supersymmetry and attempts to incorporate generational 
symmetries in the standard model (Gunion et al. 1990). We are particularly interested 
in these models because they possess stable domain wall solutions within which the Higgs 
fields may be very small. The domain wall solutions arise as an inevitable consequence of 
discrete symmetries imposed to prevent phenomenologically unacceptable flavor changing 
neutral currents. For example, if the field (pu gives mass to up quark and (po to the down 
quark, one imposes a symmetry: (pu -^ —<t>Ui Ur -^ —Ur, with Ur being the right-handed 
up quark. This prevents tree level processes in which neutral Higgs intermediate states allow 
up type quarks to change into down type quarks (Glashow & Weinberg 1977). Of course 
these symmetries mean that at least at the classical level the theories possess degenerate 
vacua, and the domain walls are stable classical field configurations interpolating between 
them. 

The domain walls are interesting to us because in their interior the gauge field masses 
may be very small and the barrier to baryon number violating processes is therefore small. 
Furthermore, as pointed out by Preskill et al. 1991, in these particular theories the domain 
walls pose no problems for cosmology, since they vacuum degeneracy is actually broken by 
very small instanton induced effects, the vacuum energy being different by ~ |/^m^ ~ (80 
MeV)*^. They may in fact be useful for producing a cosmological baryon asymmetry, as we 
discuss below. 

The electroweak baryon number violating processes violate baryon and lepton number, 
but conserve their difference. A unit change in the gauge/Higgs winding number destroys 3 
baryons and 3 leptons. Since quarks, but not leptons, are allowed to change flavors through 
Cabibbo-Kobayashi-Maskawa mixing, one can have the processes 3n ^ 17^ + T77[ + ^v, or 



2n ^n + i72 + T7^ + T\. A key prediction of our mechanism is that the high energy neutrinos 
measured from these processes would exhibit an equal excess of antineutrinos from each 
generation. 

The basic equation governing the rate of baryon number violation in the standard 

model reads, 

duB A/" , 9^ „ 9 , , , , 

— = ^(,?E.B-,?e.b). (10) 

where hb is the baryon number density, A/" = 3 the number of generations, g2 and gi the 
SU(2) and U{1)y gauge coupling constants and E, B, e, b the corresponding electric and 
magnetic fields. The linear combination occurring on the right hand side of (|I0|) vanishes 



for electromagnetic fields: baryon number violation requires the W and Z fields. If the 
wall moves at velocity f„, then the number of neutrons hitting the wall per unit time per 
unit area is 7w^AfWw, where n^ is the ambient neutron density. Requiring that a sizeable 
fraction of these neutrons are converted then leads to an estimate of the required electric 
and magnetic fields. 

Assuming equipartition between electric and magnetic fields, and a semi-relativistic 
wall velocity, we find the needed fields are: 



Brequired ~ Erequired ~ \J l^nNV^niH / a2 ~ (l.SGeV)^, (11) 

with ^2 = S'l/^TT ~ ^ the weak fine structure constant. Here we have integrated the 
anomaly across the wall width, taken to be the inverse Higgs mass itih ~ lOOGeV, 
and have ignored gi since it is much smaller than g2. The resulting fields are indeed 
large, but they are very much smaller than the electroweak scale. We have ignored the 
Cabibbo-Kobayashi-Maskawa mixing angle suppression. The B violating processes would 
presumably create heavy quark states which would be trapped upon the wall since they are 
light on it, which would annihilate with ups and downs via CKM mixing. 

Could the required fields be accreted on the wall? It is clear that the physics of the 
domain walls is quite complex and there are many possible contributing effects. Fermions 
are attracted to the wall, since their masses are lower on it, and form bound states which 
may support currents flowing without resistance on the wall. Since the neutrons carry weak 

2 

isospin, their presence automatically causes a weak isospin electric fields of order gn'fq on 
the wall if the gauge fields are massless there. But the most likely mechanism for generating 
the required ~ 1 GeV^ fields is the Meissner effect - W and Z gauge fields would be swept 
along by the back edge just because they are light on the wall and massive off it. It is also 
possible that a dynamo mechanism would operate on the wall, building up the magnetic 
fields to the saturation value -Bcrit ~ (lOOGeV)^. 
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How much energy do the fields in equation ([Tl|) cost? And how does this compare with 
the rate of energy input in the form of neutrons onto the wall? Equating the energies per 
unit area B'^equired/^^H ~ nNLjymiy where Lj^ is the path swept up in the cold neutrons, we 
find that the required field energy could be accreted in a distance as short as ~ 3 neutron 
spacings at nuclear density, or Ln ~ (a2^Ar)~^- So there is certainly no shortage of energy 
input to maintain the fields. 

We require aligned electric and magnetic fields in order to violate baryon number: we 
may think of the process occurring in two stages. First magnetic fields are established on 
the wall via some kind of dynamo mechanism coupled to the Meissner effect. In a constant 
SU{2) magnetic field, the lowest energy levels for the quark consists of left handed particles 
with positive isospin (up quarks) moving against the field and left handed particles with 
negative isospin (down quarks) moving along the field. Right handed particles do not 
couple to the field. A beam of neutrons streaming into such a field would undergo isospin 
separation, leading to an electric field aligned with the magnetic field. The resulting change 
in the gauge field winding causes neutrons to disappear. 

Some elementary steps towards establishing the viability of this mechanism would be: 
(i) finding the energy barrier required for baryon number violation in the presence of a 
domain wall, (ii) studying the fate of currents on the wall, and whether they can support 
magnetic fields of the required magnitude, and finally (iii) studying the possibility of a 
dynamo mechanism that would establish such fields out of a seed magnetic field, that is 
likely to be present in the neutron star. We shall return to these questions in future work. 

Assuming the electroweak baryon number violation processes are sustained in the 
bubble wall, the latter will be pushed forward through the neutrons by the excess pressure 
of the hot medium behind it. Presumably it will reach a terminal speed determined by 
collisions with the neutron matter. The self-sustained detonation front described in §2 
could convert all the nuclear matter mass Mtot into neutrinos and photons, releasing 
~ 5 X 10^^ ergs x (Mtot/SMo) over a time ^ (30 km/c) ~ 10"^ sec. 

When the electroweak domain wall emerges from the surface of the neutron star, the 
domain wall will still expand due to the thermal pressure behind it out to i?max/-R* ~ Pv 
where i?* is the radius of the neutron star. The lack of additional fuel will prevent it from 
expanding further; beyond this radius, the radiation fireball will diffuse out and separate 
from the domain wall. This energy release will be the primary source of energy for the 
GRB. As discussed in the next section, typically py ~ 1/100, and the limiting radius would 
be -Rmax ~ 10^ cm. If py = the wall will not expand forever, but rather at i? ~ 10^° cm 
it will contract again due to wall tension. Eventually, the wall would lose the radiation 
pressure support behind it and collapse again. The Higgs energy of the wall and the 
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false vacuum behind it will be radiated as high energy (~ 100 GeV) neutrinos, photons, 
and particle-antiparticle pairs. In contrast to supernovae, where neutrinos carry 99% of 
the energy (since photons are trapped by the opaque stellar envelope), there should be 
comparable energies in neutrinos and photons in the resulting fireball here. In addition to 
the ~ 10^^ ergs radiated in neutrinos of energy ~ 100 MeV, a second sub-millisecond burst 
carrying ~ 10^^ ergs of ~ 100 GeV neutrinos is a generic prediction of our model. 



5. Implications for Electroweak Theory and Cosmological Baryogenesis 

Electroweak domain walls can be cosmologically disastrous. When formed in a first 
order phase transition, they will come to dominate the matter density of the universe 
(Zel'dovich et al. 1975). If we assume about one domain wall per horizon volume, 
electroweak domain walls would dominate the matter density of the Universe if they 
persist until its temperature is ~ 1 keV. Thus, there is a cosmological lower bound on 
the potential difference AV = V" - 1/+ > (1 keV)^. Preskill et al. (1991) argued that 
discrete symmetries could be broken anomalously through QCD effects, yielding an energy 
difference AV ~ (80 MeV)^. The vacuum energy admits an equation of state p = —p. The 
energy density of a neutron star is ~ (200 MeV)^. Requiring the postshock pressure to be 
positive, we obtain an upper bound to this energy asymmetry AV" <^ (150 MeV)^, in order 
for the burning front to be self-sustained. This bound is well above the value predicted 
by Preskill et al. (1991). We also see that such electroweak explosions can only occur in 
matter which has densities close to nuclear, i.e. only in neutron stars. 

A related enigma in cosmology is baryogenesis. Baryon number is strongly violated 
above the electroweak phase transition, and in some theories (where baryon number minus 
lepton number is conserved) baryon asymmetries generated before this transition would 
be erased. Baryon number is conserved after the electroweak symmetry breaking, so the 
natural time to create baryon number would be during this phase transition. The Sakharov 
criteria require that baryon number B, charge conjugation symmetry C, and parity CP be 
violated, and that the system be out of equilibrium. In the minimal standard model of 
particle physics, two of these criteria are problematic. The current measurement of the top 
mass rritop ~ 180 GeV appear to preclude a first order phase transition, which would be the 
natural way to bring the Universe out of equilibrium. Also, the magnitude of CP violation 
in the minimal standard model cannot produce the observed baryon-to-photon ratio 
f] ~ 10"^*^. Moving topological defects such as domain walls restore the out-of-equilibrium 
requirement (Prokopec et al. 1996), and the extra freedom in selecting the second Higgs 
field also allows a match for the required CP violation rates. The long persistence of 
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the anomalous domain walls will also amplify the net baryon production to become more 
consistent with the small observed CP violation rates. 

Another consequence of the long persistence of these walls is a coherence length of the 
baryon production regions which is intermediate to the horizon size at electroweak and the 
horizon size when the domain walls disappear. This length scale of ~ 10^^-10^^ comoving 
cm, is significantly longer than the diffusion length at nucleosynthesis (Jedamzik and Fuller 
1995). The baryon-to-photon ratio would generically fluctuate in different regions, resulting 
in inhomogeneous nucleosynthesis and affecting the standard predictions for the abundance 
of light elements. 



6. Conclusions 

We have argued that the gravitationally unstable core of a neutron star merger might 
trigger an electroweak burning front which would propagate outwards and convert all the 
nuclear material into radiation. The associated detonation front would be self-sustained and 



propagate at a constant semi-relativistic speed ~ J?>/Ac). The baryon-number violating 
reaction is induced by a domain wall, which is naturally predicted in double Higgs models 
of electroweak interactions. We have shown that the electroweak burning front can only 
be self-sustained if it propagates into material with nuclear densities. Although our 
mechanism provides an efficient way for alleviating the baryonic contamination problem, 
some surrounding stellar debris might remain unburnt and limit the Lorentz Factor of the 
expanding fireball. Minimal contamination is required to make the emission spectrum from 
the fireball non-thermal (Goodman 1986; Paczyhski 1986). 

The electroweak explosion process releases ;:^ 10^^ ergs over ^ msec, and produces 
a fireball which subsequently impacts on the ambient (interstellar) medium. This 
impact generates an external shock which could produce via synchrotron emission and 
inverse- Compton scattering of Fermi-accelerated electrons, the observed spectrum of GRBs 
(Meszaros & Rees 1993a,b; Meszaros, Laguna, & Rees 1993). In this model, the observed 
variability of GRBs might be either due to inhomogeneities in the ambient gas density, or 
due to instabilities at the shock front (e.g. Waxman & Piran 1994). The release of ^ 10^"^ 
ergs remedies the energy problem which was previously identified for external shock models 
(Sari & Piran 1997a,b), and is in fact required if GRBs follow the star formation history 
of the Universe with the dimmest bursts originating at redshifts as high as 2; ~ 6 (Wijers 
et al. 1997). The integrated light from afterglow observations can set a lower bound 
on the total energy released, and in principle challenge models which provide too little 
energy. The release of ~ 10^^ ergs in more conventional astrophysical settings (see e.g., the 
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hypernovae model of Paczynski 1997) is likely to be contaminated by baryons at a level 
which depends on many random parameters, such as the orientation of line-of-sight through 
the source. Under such circumstances, one would expect a continuous distribution of events 
with varying levels of baryonic contamination, from the level of ^ 10~'^Mq required for 
GRBs and up to the level of ~ IOMq found in supernovae. Such models would therefore 
predict the existence of softer events, such as X-ray bursts (which were not observed so 
far), in addition to GRBs. Our model naturally accounts for the gap between the baryonic 
contamination levels found in GRBs and supernovae. 

The collapse of the bubble wall could also release a fraction py of the neutron star 
rest mass, i.e. ~ 10^^ ergs, when some of the Higgs potential energy is converted into 
non-thermal neutrinos with very high energies ~ 100 GeV. The number of ~ 100 GeV 
neutrinos per unit area from a source at a distance ~ 3 Gpc would be ~ 10® km~^. Full 
timing and directional coincidence with GRB events might allow statistical detection by 
large area Cherenkov arrays, such as AMANDA. With a typical conversion efficiency of 
10~^ into upward moving muons (Gaisser et al. 1995), a year long cross-correlation with 
~ 300 bursts could yield ~ 3 events with a km^ array. This neutrino burst supplements the 
predicted neutrino flux from standard ffieball models (Waxman & Bahcall 1997, Paczynski 
and Xu 1994). A larger amount of energy ~ 10^^ is released in ~ 200 MeV neutrinos, 
remnant of the burnt phase behind the detonation front; however, these lower energy 
neutrinos would be more difficult to detect. 

Finally, we note that the extended Higgs models are amongst the simplest extensions 
of the standard model, and the discrete symmetries needed to produce the domain walls are 
required by particle physics phenomenology. These models may well provide a new solution 
to the problem of electroweak baryogenesis, with associated inhomogeneous cosmological 
nucleosynthesis. The quantitative constraints which are placed on the double Higgs theory 
by our GRB model could be tested with future accelerator experiments. Conversely, if our 
proposed link between electroweak theory and GRB is established, GRBs might be used 
to place constraints on parameters of the Higgs sector which go beyond the capabilities of 
laboratory accelerators. 

We thank Shmulik Balberg and Eli Waxman for useful discussions. This work was 
supported in part by the NASA ATP grant NAG5-3085 and the Harvard Milton fund (for 
AL and UP). 
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Fig. 1. — Self-Similar explosion solution for a = and pv = 0. The different curves 
correspond to different detonation front speeds, A^. The solid line is the generic solution in 
the presence of weak viscosity. 



